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Abstract. The presence of impurities is considered in gyrokinetic calculations of 
ion temperature gradient (ITG) instabilities and turbulence in the reversed field 
pinch device RFX-mod. This device usually exhibits hollow Carbon/Oxygen profiles, 
peaked in the outer core region. We describe the role of the impurities in ITG 
mode destabilization, and analyze whether ITG turbulence is compatible with their 
experimental gradients. 
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On the mutual effect of ion temperature gradient instabilities and impurity peaking in the reversed field pinch'. 
1. Introduction 

Impurity transport in magnetic fusion devices is a long-time debated subject, which in 
diverted tokamaks gains topical relevance to avoid high effective charge accumulation in 
the plasma core. This has been a rich field of research in experimental fusion science for 
many years, with remarkable attention to the role of auxiliary heating systems on the 
impurity flux (often reversing impurity convection from inwards to outwards) [H [21 [3l H] , 
and to the role of edge [5] and internal [6] transport barriers in inducing impurity 
flattening or accumulation. Much of the related theoretical activity has concentrated 
on the role of neoclassical transport [7], and more recently on that of electrostatic 
ITG/TEM turbulence [H HI [TOl [H], but often with contrasting results. Viceversa, 
whether large impurity concentration may yield turbulence enhancement or quenching, 
is an even more open issue. In the context of electrostatic turbulence the reversed field 
pinch (RFP) is a scarcely explored configuration - transport being historically thought 
as dominated by large-scale MHD turbulence - which may give a contribution to the 
knowledge of impurity transport in toroidal fusion plasmas. 

Recently the study of microturbulence has in fact extended to the RFP 
configuration. The main reason is the achievement of core plasmas naturally endowed 
with almost conserved magnetic fiux surfaces. In the RFX-mod device [12j, for example, 
this has been accomplished thanks to the self-organized single helicity states, in which 
a single m = 1 kink-tearing mode is able to dynamically sustain the discharge. The 
net effect is the arising of an overall helical warping, usually with the periodicity of 
the innermost m = 1 resonant mode. When the plasma current is high enough, the 
plasma has a further transition towards the so called single helical axis state [13], with 
separatrix expelled and a rather peculiar electron temperature profile, fiat in the core 
and barrier-like at mid radius. Present diagnostics at RFX allow for determining the 
electron temperature only, but measurements from passive spectroscopy let estimate for 
the electron to ion temperature ratio Tg/Tj a value about 3/2 in the core. 

Does microturbulence really play an important role in RFP transport? Previous 
studies have shown that ion temperature gradient (ITG) instabilities [HI [15] can be 
considered only a partial contributor to particle/energy transport in the RFP [T6l \T7\ 
[T8l [To] , microtearing modes being much more relevant [20]. This has been explained 
as a consequence of the significant Landau damping of the ITG mode, driven by the 
field connection length ~ 27ra (a is the torus minor radius), or equivalently by the 
low value of the safety factor, g ^ 1. Part of this paper is devoted to a more careful 
analysis of the ITG stability threshold. In particular we will investigate the role of 
the impurities in ITG mode destabilization. Such a destabilization may occur in RFX- 
mod, due to the hollow profile of the impurities [21] caused by their inefficiency in 
reaching the inner plasma [22]. For this analysis we adopt the code GS2 [231 121], 
a widely known gyrokinetic code based on the electromagnetic nonlinear gyrokinetic 
equation [211 [211 [22], solved in a fiux tube domain. The modifications required for GS2 
to work in the RFP configuration are examined in Ref. [17], and mainly concern the 
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large poloidal component of the magnetic field. The RFX-mod single helicity fields have 
been approximated as axisymmetric: 3D gyrokinetic simulations based on reversed field 
pinch VMEC equilibria [2H1 EH] is an upgrade to look at for future studies. 

In RFX-mod the most populated impurity density profiles are hollow, peaked in 
the outer core domain 0.6 ^r/a ^ 0.9. This feature is encountered in all the RFX-mod 
experimental regimes, either during the chaotic multiple-helicity stages or in the laminar 
quasi single helicity states [22j. The absence of stationary impurity accumulation is an 
appealing feature for any fusion device, which makes this RFP machine an interesting 
test bed both for experimentalists and theoreticians. Since neoclassical transport is 
practically negligible in this configuration, impurities turn out to be driven either 
by some magnetic chaos in the outer core region or, possibly, by small scale micro- 
turbulence, an option never explored so far. In this work we are going to look into the 
latter possibility, neglecting the effects related to magnetic surface destruction due to 
tearing or microtearing instabilities. In particular we will consider the role played by 
purely electrostatic ITG turbulence. Dealing only with this kind of turbulence can be 
heuristically considered relevant if the resulting ion thermal diffusivity is much larger 
than the cross-field ion diffusivity due to magnetic chaos. Assuming a gyro-Bohm scaling 
for the former (~ pfvih^i/Lxi, with pi ion Larmor radius, Vth,i ion thermal speed, Ljv 
ion temperature gradient scale length) and a collisionless quasi-linear expression for 
the latter (~ b'^LcVth,i, with b normalized magnetic field perturbation and Lc parallel 
correlation length), such condition roughly reads {pi/a)'^{a/ LxJ/b'^ ^ 1, where we have 
considered Lc ~ a, a valid estimate for both tearing and microtearing turbulence in 
RFX-mod. An alternative to get the same condition is to require the ITG mode 
growth rate to be much larger than the rate needed for the ions to spread radially 
by ~ 1 — 10 Pi (~ ^x\tg) magnetic diffusion. The inequality above is easily 

fulfilled in the case of micro-tearing turbulence, where b ~ pe/Lx^', conversely, for a 
large scale tearing turbulence b is much larger, often yielding {pi/ aY{a/ Lt^) /b'^ < 1. 
Thus, as regards ion transport, neglecting magnetic macro-turbulence turns out to be 
the strongest assumption. Indeed, as we will see in the final part of the paper, the action 
of non-ITG transport mechanisms needs to be invoked to explain the observed impurity 
profiles. 

In the following we describe the effect of impurities on ITG mode destabilization 
in Sec. [21 then the role of ITG instabilities in determining impurity transport in Sec. |3l 
and finally the nonlinear mutual interaction between them in Sec. HI 

2. Effect of impurities on ITG modes 

In RFX-mod the density profiles of the most populated impurity ions are usually 
peaked outside the core, as shown in Fig. [H Carbon and Oxygen being the most 
relevant species. Based on passive spectroscopy evaluations, their volume content is 
evaluated to be around a 1-2% of the total number of electrons. The present Carbon 
and Oxygen profiles are reconstructed following the method described in Ref . [30] . All of 



On the mutual effect of ion temperature gradient instabilities and impurity peaking in the reversed field pinch 

the available experimental data concerning emission lines, continuum radiation and total 
radiated power emitted by the intrinsic impurities are reproduced with a one-dimensional 
(radial) time- dependent collisional radiative diffusion model. The model gives the best 
evaluation of the impurity transport coefficients with the resulting Carbon and Oxygen 
ion density profiles. Such profiles are fully consistent with the output of the transport 
code RiTM [31] or its MHD extension Riport |32]. Given the electron and impurity 
density profiles, the ion density is derived accordingly. 

In the following of this section we are going to describe the setup of the gyrokinetic 
simulations, the effect of a single impurity species on ITG modes, and finally the effect 
of multiple impurity species (CV-Vll and OVll-IX) on these modes. The reference 
experimental cases are the well-examined shots 23977 at 200 ms and 24932 at 64 ms. 
Both are single helical axis cases which experience a barrier in the electron temperature 
profile at mid-radius. The first one has the electron density peaked in the core, while 
the second has a slightly hollow He profile until p = r/a~2/3. In shot 23977 (cf. Fig.[T]) 
the most abundant Carbon ions in the outer core, CV (or C4+) and CVI (C5+), have 
logarithmic derivatives ajLn^^ = —a{dncy / dr) / ncy — —12 and a/Ln^^^ ~ —6 in the 
radial interval 0.6 < r/a < 0.7; similarly the most abundant Oxygen ions, OVII (or 
06+) and OVIII (07+), have ajL^Q^^^ — —11 and a/L^Q^^^^ ~ —5 in the same region. 
In shot 24932 the impurity radial profiles are even steeper, their logarithmic gradients 
locally reaching values ^ —15 for CV and OVII, in r/a G [0.6,0.9]. 

2.1. Background 

Before discussing the setup and the results of the simulations, we briefly review the 
modifications required for GS2 to work in reversed field pinch geometry. We recall 
that such geometry is characterized by a g profile everywhere -C 1, decreasing from the 
core to the edge, where it becomes slightly negative. As extensively discussed in [T7] . 
having the poloidal and toroidal component of the magnetic field the same strength. 
Be ~ the definition of some terms in the gyrokinetic equation has been modified 
with respect to the standard tokamak geometry: (i) the curvature and V-B drifts in the 
drift frequency Ud include terms of order of 0{1) proportional to Bl/ p and —OB /dp 
respectively; (ii) the parallel gradient is b ■ V oc {Bq/ p)de, where the poloidal angle 9 
is used as longitudinal (ballooning) coordinate along the flux tube; (iii) the equilibrium 
field is given hy B = Bq{1 — 5 cos 9) with 6 < ep, e = a/R; (iv) the binormal component 
of the wavenumber is given by ky = {n/r){e^p'^ + q'^Y^'^, where n is the toroidal mode 
number. 

Coming to the description of the quantities related to particle distributions, 
we recall that an arbitrary number of species can be included in GS2, with a full 
local description of charge zj, mass mj, density uj and density logarithmic gradient 
a/Ln- = —a{dnj/dr)/nj, temperature Tj and temperature logarithmic gradient a/ Lt^ = 
—a{dTj/dr)/Tj, and coUisionality uj. Plasma quasi-neutrality has to be preserved across 
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In this work most of the hnear simulations are performed including a single 
impurity only, with charge z and density n^, so that the expressions above simply read 
Ue = Ui + zuz, and ajL^^ = {ni/rigja/ + {znz/ne)a/ Ln^. The electron density and 
temperature are the experimental ones. Consistently with the RFX-mod experimental 
operation, a H background is considered, = 1. The conditions Te/T, = 1.4 and 
Tz = Ti are assumed Vp G [0, 1]. Simulations are performed collisionless, which reduces 
the computational time, since the effect of electron collisions has been tested on a few 
representative cases and found to be small. Trapped electrons are always included. 
Concerning the impurity density profiles, these have to be considered as "starting" 
profiles: in the following we will in fact scan over gradient and density, so as to provide 
a sensitivity study beyond the experimental uncertainties. 

The way we carry out single impurity calculations is to fix one impurity ion 
at a time (among the most relevant ionization states), and to rescale its density 
profile so that the volume concentration is a 2% of the total number of electrons, 
Jq Uzrdr/ Jq Uer dr = 2 ■ 10~^. Of course this implies an overestimate of the density of 
the ion under consideration, but allows to have lighter gyrokinetic simulations and to 
better understand the physical mechanism at work. 

As a first result, in Fig. |2] we show the ITG stability study taking into account 
CV (the most populated state of Carbon) or CVll (fully stripped Carbon), for shot 
23977. The stability threshold is calculated as the minimum value of the ion temperature 
gradient at which the growth rate of the mode is positive. The comparison with the case 
without impurities shows that the inclusion of an impurity species with positive radial 
derivative (a/L„^ < 0) causes ITG mode destabilization. Conversely, a negative radial 
derivative (a/L„^ > 0) induces ITG mode stabilization. A higher stability threshold 
implies an overall lower growth rate of the ITG modes, and viceversa, while the real 
frequency remains almost untouched, see Fig. [2]-c. For the present profiles the effect of 
the inclusion of a single impurity is not dramatic: in either cases the thresholds remain 
much larger (a factor ~ = R/a) than for a tokamak. Landau damping thus remains 
the key mechanism discriminating between the two configurations. As we will see, a 
larger impurity density derivative is required in order to get an important ITG mode 
excitation, or even the occurrence of impurity drift instabilities. 

2.2. Parametric study with a single impurity species 

To understand the interplay between impurities and ITG modes, we extend the single 
impurity analysis to a larger parameter domain. In particular we focus on the 
dependence of the mode frequency /growth rate and stability threshold on the local 
quantities nz{p)/ne{p) and a/LnXp)- Based on the experimental case of Fig. [2l we 
perform a set of linear calculations at p = 0.7 for CV, and at p = 0.6 for CVII. The 
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local impurity density is varied in the interval Uz/ue € [0, 3%], and the gradient in the 
interval a/L„^ G [-20, +20]. 

The result of this set of computations is reported in Fig. [3l where the critical ion 
temperature gradient is shown as a function of both the parameters of interest. Large 
impurity densities and positive radial gradients of impurity density profiles (a/L„^ < 0) 
are favorable to ITG destabilization. If the whole impurity content at p = 0.7 were 
CV (Fig. Ha, "CV @ 2% vol. concentr.") or CVII at p = 0.6 (Fig. [3]-b, "CVII @ 2% 
vol. concentr."), the destabilization would not be so relevant. However, in the case of 
CVII, we can see that for high enough n^/ng and a/L„^ ^ —1, ITG modes are well 
destabilized, and mode destabilization can be found also for a vanishing ion temperature 
gradient. The effect is visible only in the CVII case, due to the higher charge of the ion 
and to the consequent higher charge fraction zriz/rie- The resulting instability is the 
so-called impurity mode, well discussed in the literature [33l [HI EHl [36l [371 US]- The 
effect of the impurity profile becomes even larger for larger ion charges, e.g., for the 
most ionized states of Oxygen. 

Let us investigate more quantitatively the effect of an increasing local impurity 
concentration at fixed outwardly peaked impurity, a/L„^ = —20, for the case CVII at 
p = 0.6, Fig. [H Increasing Uz/ue leads to a progressive destabilization of the ITG 
mode, until the growth rate 'y{a/LTi) becomes positive for a vanishing ion temperature 
gradient. In the inset of Fig. [Hthe resulting impurity mode structure at Uz/ue = 3% as a 
function of the ballooning angle does not reveal a major modification with respect to the 
usual ITG mode, neither the mode frequency has discontinuities at low ion temperature 
gradients. As reported in Ref. [37] the impurity mode and the r]i mode (i.e., the ITG 
mode, rjj = LnjLT^ henceforth) are strongly coupled, making it difficult to distinguish 
between them. 

As a general comment on the effect of the densities on ITG instabilities, for a 
fixed Ue profile the shape of the impurity density has strict consequences on the main 
ion density profile and derivative. We recall that in the RFP the (hydrogenic) ITG 
threshold a/IvT,, critical as a function of a/Lm has a minimum for a/L„. > [a/Lm — 4 
in this case), mainly due to the strong cjvb resonance occurring around such density 
gradient, cf. Fig. 6-b of Ref. [l7]. Therefore, given the experimental Ue profile, in order 
to preserve quasi-neutrality, a negative a/ Ln^ (outwardly peaked impurity) causes a/L„. 
to increase, and such increase (unless very large) provides ITG mode destabilization. 
Viceversa, positive a/L„^ (inwardly peaked impurity) causes a/L„. to decrease, and the 
ITG mode to be progressively stabilized. 

Looking again at Fig. [HI for the experimental almost flat Uf, profile (|a/L„^| < 1), 
if the impurity density is high enough, the effect of a large impurity density gradient 
is to cause an inversion of the sign of the main ion density gradient. In both cases, 
CV and CVII, the (left) half-plane a/L„^ < roughly corresponds to a/Iv„. > 0, 
with a maximum for a/L„. corresponding to the upper-left corner of the respective 
frames {a/Ln^ ~ 3 for CV, a/Ln^ ~ 4 for CVII). Conversely, in the (right) half-plane 
a/Ln^ > the value of a/L„. becomes negative {rji < 0), with the minimum in the 
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upper right-corner. This study has been repeated for other ionized impurities, shots and 
electron density profiles, indicating that for a strong ITG mode destabilization or even 
the impurity mode onset to occur, there must be a combination of the destabilizing term 
a/L„. > and the relative sign of impurity and main ion density gradients, oJ^iU^z < 
(with oj^j = kycTj / eB Lnj) , cf. Ref. [33]. In particular the impurity mode ai rji = = rj^ 
is produced for reversed impurity density gradients, and if the impurity charge fraction 
is high enough. This phenomenology confirms the former results obtained in tokamak 
plasmas [331 EH ESI ESI EZ], the physical mechanisms at work being essentially the same 
in both cases: however, the manifest drop of the stability threshold observed in the RFP 
makes impurities particularly important, as they can determine the transition from a 
condition where ITG modes are inherently absent to a scenario possibly dominated by 
such turbulence. 

The parameter region providing impurity-mode destabilization is not reached in the 
experimental cases examined so far. Such results are obtained by modifying the single 
impurity parameters of interest. To get a more realistic conclusion we have to consider 
the case of many partially-ionized impurity ions simultaneously present at a time, as in 
the experiment. 

2.3. Simulations with several impurity species 

The 6 most ionized ions. Carbon V to VII and Oxygen VII to IX, are included in 
the linear ITG simulations, overall rescaled so as to get a total 2% concentration, 
e lo ^''^f jQ^e^dr = 2 ■ 10^^. Again, for any impurity species the colhsionality 
is neglected, Uj = 0, and the respective temperature profile is assumed equal to the ion 
one, Tj = Ti, a/Lx^ = a/Lx^. The quasi- neutrality across the flux tube is imposed to 
the main ion density, making Eq|T] satisfied all along the radius. 

We would expect a cumulative effect of such impurities if their profiles had all either 
a positive or a negative slope, but this is not generally the case. The single impurity 
profiles are not peaked at the same radius, and the stripped CVII and OIX profiles are 
actually rather fiat. This is for example the case for shot 23977, and also for the other 
shot we are going to analyze now, 24932 at 64 ms. 

In Fig. |5] the two shots are compared. They are characterized by different electron 
density profiles, the latter being slightly hollow. The net effect of the impurities is in 
both cases a somewhat "destructive interference" . No dramatic ITG destabilization is 
found in the region of interest, neither in the peaked (23977) nor in the hollow (24932) 
electron density profile case; only a slight destabilization is present at p ~ 0.7, where 
a/^n, > 1. 

We have performed another set of linear gyrokinetic simulations based on a different 
electron density profile. Starting from the data of shot 23977, and in particular keeping 
the same impurity profiles, we have built a peaked density, as shown in Fig. [6]-a. 
Impurity profiles are rescaled on the new rig profile to get a total 2% concentration. 
A similar electron density is not strictly experimental, peaked profiles emerging only in 
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very low density discharges, but remains a useful test-bed for our purpose. Following 
exactly the same procedure carried out before, we look at the different characteristics 
of ITG modes without impurities and with 6 impurities (CV— i-VII and OVII— i-IX). 
While the introduction of impurities does not provide any significant changes in the 
core until p ~ 0.6 (Fig. EJ-b/c), in this case their superposition causes an important 
mode destabilization at p = 0.7, where a/L^^ ~ 4 and C(;^,jC(;*2 < for the most abundant 
impurities. Such a destabilization is obtained also for vanishing rji, see Fig. [6]-d. The 
spectrum of the resulting impurity mode at a/Lx^ = is plotted in frame f, and covers a 
larger kyPi-iange with respect to the usual ITG mode excited, for example, at a/Lx^ = 6. 
At the extreme of the radial interval, p = 0.8 (Fig. |6]-e), the dependence of the growth 
rate as a function of a/Lx^ appears less simple: the values of 7 for low a/Lx^ are 
missing because we are not considering here instabilities with negative real frequency, 
i.e., propagating in the direction. In the corresponding range of a/Lxi (covered by 
the shaded area in the figure), trapped electron modes appear to be the fastest growing 
instabilities, indeed excited by the high value of the electron density gradient at that 
radius. The occurrence of trapped electron modes in the reversed field pinch is a topic 
to be extensively investigated in the future. 

Summarizing this first part of the paper, no dramatic influence of the impurities is 
found on ITG mode destabilization for rather flat electron density profiles. In the next 
section we will discuss a method to indirectly evaluate the nature of the microturbulence 
at work, based on the direction of the impurity flux. This could be considered as a 
counter-evidence of the impurity- driven nature of microinstabilities. 



3. Effect of ITG modes on impurities 

A completely different approach is used to evaluate linearly the direction of the 
impurity flux in ITG dominated discharges. While for the previous case the impurity 
concentration had to be large enough to substantially perturb the system, for this study 
impurities are assumed to have negligible charge fraction, zn^/ne <^ 1, so as not to 
influence the turbulence itself. In other words, we are supposed to deal with trace 
impurities. 



3.1. Background 

We summarize how, from linear gyrokinetic studies, it is possible to deduce the so-called 
peaking factor of an impurity species. The procedure is more extensively described, for 
example, in Refs. [8], [38], and [39] . 

For any species j, the quasi-linear particle flux Tj = {fijVExB), with (■) magnetic 
flux surface average, can be written in the form 
oF, ^ a ^ a 

-^ = D,— + Dtj — + avpj (2) 

Uj Lit J 

with Dj diagonal diffusion coefficient, Dtj off-diagonal thermodiffusion coefficient, and 
vpj pure convection. In general, these coefficients do depend on ajL^^ and ajL^^ through 
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the dependence of the eigenvalues on such gradients. In the case of a single impurity 
with charge 2;, if its charge concentration is so small {zUz/Uf. ^ 1) that it can be 
neglected in the gyrokinetic Poisson/ Ampere equations, the coefficients Dz, Dtz and 
vpz become independent of a/L^^ and ajLrp^. Thus the flux of a trace impurity species 
simply becomes a linear function of the gradients. 

The equation above for the impurity flux can be written as 



DA-^^ Ct-t- + Cp (3) 

where we have defined Ct = Dtz/Dz and Cp = avpz/Dz. Under the mentioned 
assumption of trace concentration, such coefficients can be derived from the slopes of 
the quasi-linear flux Tz (calculated by GS2) as a function of the logarithmic gradients 
a/Ln^ and a/Lx^, and from the residual flux at vanishing gradients. The way to get 
them is to turn on/off the gradients, always preserving quasi- neutrality. The total 
convective contribution is flnally given by the sum of the thermodiffusive term and the 
pure convection, 

^ = Ct-^ + Cp. (4) 

In stationary conditions, and without particle sources, diffusion and convection must 
balance so as to yield = 0. This reflects on a constraint for the impurity density 
gradient, a/L^^ = —Ct a/ Lt^ ~Cp = —aVz/Dz- The sign of the coefficients Ct and Cp 
and that of the peaking factor —aVz/Dz are strongly affected by the dominant instability, 
and can be considered the signature of the related turbulence in the plasma. 



3.2. Impurity peaking in ITG dominated plasmas 

Following the procedure just described, we now study how the impurity flux is affected 
by ITG instabilities. In the absence of plasma rotation ITG modes are commonly 
believed to yield an inward impurity flux, at least in tokamaks [U HQ], HH [11]. The RFP 
turns out not to be an exception to this rule. 

For our set of linear simulations the charge density is set to zuz/n^ = 10~^, so as 
not to signiflcantly perturb the system. The ITG instabilities are artiflcially excited by 
imposing an ion temperature gradient above the threshold at every radius in the region 
of interest, in this case a/LT^ = 2a/ Lt^. This is not consistent with any proflle, but 
is certainly allowed in the framework of flux tube simulations. 

The dependence of Ct, Cp, and aVz/Dz on the trace impurity charge Z, at flxed 
radius, is shown in Fig.O-a. The impurity charge is varied between 6 and 28, i.e., between 
Carbon and Nickel. The frequency is positive and constant. The mode structure itself 
is not changed, indicating that the ITG mode is not modifled by the inclusion of a 
low-density impurity. The t her mo diffusion is positive and decreases with Z, with an 
opposite sign with respect to the curvature pinch Cp. The quantity aVz/Dz turns out to 
be negative, i.e. an inward pinch is driven; analogously, the peaking factor is positive, 
resulting in an impurity accumulation in the core and a peaked impurity proflle in 
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stationary source-free conditions. This is not compatible with the experimental results 
discussed in Ref. [22], which report a strong outward impurity convection, especially in 
the outer plasma core. 

The impurity convective contribution aVz/D^ is also evaluated across the ITB, at 
several radii for a fixed charge, e.g.. Nickel, see Fig. [3-b. Until the dominant instability 
is ITG, the inward direction of the pinch cannot be prevented. Of course, the resulting 
radial trend of aVz/D^ has to be merely seen as a qualitative indication, being strictly 
dependent on the selected local ion temperature gradient, which is actually unknown; 
however, the sign of the peaking factor is well-established. Assuming microturbulence 
to be the dominant transport mechanism, other instabilities are likely to be at work. 

4. Nonlinear studies of ITG turbulence with impurities 

Although linear gyrokinetic calculations are essential to determine the impact of 
impurities on ITG mode destabilization, and helpful to derive quasi-linearly the basic 
particle transport properties of a trace species, the problem of the mutual interaction 
between ITG modes and impurities should be rigorously treated in a self-consistent 
nonlinear framework. In the final part of this paper we start addressing the topic of 
ITG turbulence in 3-species RFP plasmas, reconsidering two cases discussed in Sect. [3] 
and [2] respectively. 

Previously we have claimed that trace-impurity fluxes are linear in den- 
sity/temperature gradients. As a numerical proof, we have performed a set of nonlinear 
simulations for Nickel at p = 0.5, with exactly the same parameters used for Fig. [71 
Only the impurity density gradient is varied in this case. Thermoconductivity is there- 
fore included in the convective term, with the impurity temperature gradient set to 
a/Lx^ = a/LTi — 8.3. At p = 0.5 the magnetic shear is s = —0.64. The wavenumber 
grid has Ak^Pi — AkyPi = 0.1, max{kyPi) = 0.6, with a box size Ly ~ 64 pj, L^/Ly ~ 1. 
The peak of the linear spectrum is at kyPi ~ 0.3. The result of the nonlinear calculations 
is shown in Fig. [HI where the perfect linearity of the flux in a/L„^ yields aVz/D^ — —3, 
i.e. a larger value (in absolute value) than the linear results obtained in the previous 
section, but not very far from them. 

It is widely more interesting to come back to the single- impurity simulations with 
large Uz/ue of Sec. [21 as an integration to the previous linear analysis. To this aim, we 
deal again with shot 23977, Carbon VII at p = 0.6 with a/L^^ = 10, and in particular 
for Uz/ue = 0.03 with impurity gradients in the interval a/Ln^ G [—10, +10], see Fig. [31-b 
for a comparison with the linear cases. The electron density gradient is a/Ln^ = 0.23. 

At Uz/ue = 0.03 impurities do perturb the system, the ITG stability threshold, 
and the possible onset of impurity drift modes. Linear ky spectra are different in the 
considered range of a/L„^, being more extended for a/Ln^ = —10 (until kyPi = 0.7) 
than for a/Ln^ = 10 (until kyPi = 0.5). For the nonlinear simulations the wavenumber 
grid has Ak^Pi — Akypi = 0.075, with ma.x{kyPi) ~ 1 in the case a/L^^ = —10, and 
with a box size Ly ~ 80 pi, L^/Ly ^ 1. At this radius the magnetic shear is s = —1.22. 
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In Fig. [U]-a we report the temporal evolution of the ion conductivity Xi = 
—Qi/iP'idTi/dr), where Qi is the ion heat flux, in gyro-Bohm units, xgb = pf^th.i/*^- 
We clearly see the expected decrease of Xi fo^' increasing a/Ln^, in analogy with the 
increase of the stabihty threshold. In Fig. [9]-b the saturated particle fluxes are reported 
as a function of a/L„^. We notice that: (i) outwardly peaked impurities (a/L„^ < 0), 
with inwardly peaked main ions due to quasi-neutrality (a/L„. > 0), are found to 
have a net inward flux (F^ < 0), and to provide an outward ion flux (Fj > 0); (ii) 
for fiat impurities (a/L„^ = 0) the impurity flux is inward (F^ ^ —0.06) as well 
as the main ion flux (Fj < 0); (iii) inwardly peaked impurities (a/L„^ > 0), with 
outwardly peaked main ions due to quasi-neutrality (a/L„^ < 0), have an outward flux 
(F^ ~ 0.15), and cause ions to be pinched (Fj < 0). In all the cases electrons are 
pinched, Fe < 0. The anomalously large value [F^l obtained for the reversed Carbon 
gradient (a/L„^ = —10) can be interpreted as the result of the impurity mode branch 
working in the high wavenumber range, which indeed has the consequence of producing 
further inward convection of Carbon and an outward flux of Hydrogen |34] . 

As for the previous linear simulations, this rather strong effect is associated with 
the high impurity charge concentration and to the fact that the whole impurity content 
is moved onto a single ionized species, C VII in this case, with steep density profile. 

5. Conclusions 

The presence of the impurities has been considered in the framework of linear and 
nonlinear gyrokinetic simulations, applied to the RFP configuration. ITG modes, whose 
stability threshold a/Lj'^ critical in a two-component plasma is a factor ~ e^^ larger 
than in the tokamak, turn out to be slightly destabilized in the presence of outwardly 
peaked impurity profiles. This behaviour is in agreement with the past tokamak and 
RFP results. Parametric studies show that a more serious destabilization, including 
the appearance of impurity drift modes, could be obtained for very peaked main gas 
densities, a condition not commonly observed in correspondence to the positive impurity 
slopes of the outer core plasma. 

On the other hand, if ITG instabilities/turbulence only were active, the convection 
of impurities should occur towards the plasma core. Instead, experimentally it is found 
to be outward. This suggests that other instabilities may be active; further work is 
required to understand their existence and their role on impurity flux reversal. In 
particular the presence of magnetic-surface breaking instabilities has to be considered; 
even in single helicity states magnetic field stochasticity could drive a significant fraction 
of the impurity transport in the region investigated in this paper, as well as of the main 
gas particle and energy transport. 
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Figure 1. Impurity profiles of Carbon (left) and Oxygen (right), shot 23977 at 200 
ms. The total volume impurity concentration is about 2% of the total electron content. 
The occurrence of ITG modes will be analyzed in the shaded radial region. 




Figure 2. ITG threshold study for shot 23977. In (a) and (b) only CV and CVII 
are considered respectively, their volume concentration being set to a 2% of the total 
electron content; ♦ symbols represent the thresholds calculated without any impurities. 
In (c) real frequency (solid) and growth rate (dashed) spectra refer to the most external 
radius at the highest ion temperature gradient of the scan. 
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Figure 3. a/^Ti, critical contour as a function of the local values of nz/rif, and a/L„^ 
for a single impurity at fixed radius, shot 23977. At p = 0.7 (a) the experimental 
value of CV is increased to get it covering the whole 2% of the impurity volume 
concentration, then the analysis is extended to (rij/ne, a/L„^) £ [0,0.03] x [—20,20]. 
At p = 0.6 (b) the experimental value of CVII is increased to get it covering the 
whole 2% of the impurity volume concentration, then the analysis is extended to 
{uz/ne, a/ LnJ 6 [0,0.03] x [—20,20]. In the green area the instability is obtained 

with a/ Lt,, critical < 0. 




Figure 4. Section a/Ln^ = —20 of the contour in Fig. [3]-b. Increasing rizln^ drives 
destabilization also for vanishing a/Lxi- In the inset the mode structure is plotted as 
a function of the ballooning angle at a/Lxi — 0. 




Figure 5. Electron density for shots 23977 and 24932, and respective radial 
dependence of the ITG threshold without impurities (♦ symbols), and in the presence of 
the 6 most ionized impurities, CV^VII and OVII— >IX (■ symbols). The experimental 
profile of a/L^^ is plotted for comparison. 
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Figure 6. Modified peaked electron profile for shot 23977 (a), with the maximum 
growth rate 7 as a function of a/Lx. from p = 0.5 to p = 0.8 (b-e) for the case 
without impurities (♦ symbols), and with the most stripped impurities, CV— >VII and 
OVII^-IX (■ symbols). The shaded area in (e) represents the range in o/Lt; where 
modes do not propagate in the ion diamagnetic direction. Two spectra at p = 0.7 
are in frame f, for a/Lj^. = 6 without impurities (♦) and for a/ L^i — with all the 
impurities (■) 
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Figure 7. Impurity peaking av/D with respective Cp, Ct and ITG mode frequency 
Lo for different values of impurity cfiarge at fixed radial position p ~ 0.5 (a), and 
for different radial positions at fixed charge Z = 2d> (b). The unphysical condition 
Lt^cxp/ = 2 is imposed to the ion temperature gradient at each radius. 




Figure 8. Nickel flux as a function of its density gradient (with zriz/n^ ~ 10~ , 
z = 28) at p = 0.5, from three-species nonlinear simulations of ITG turbulence 
assuming Ly^oxp/iTi — 2. The impurity flux is linear in a/L^^, and av^jD^ — —3. 
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Figure 9. Three-species nonlinear simulations of ITG turbulence with CVII having 
nz/ne = 0.03, a/Ln, G [—10,10] (cf. Fig. |3]-b), and with ion temperature gradient 
a/Lxi — 10. The ITG destabilization for negative a/Ln^ corresponds to an increase 
of the ion thermal conductivity (a). Outwardly peaked impurity profiles (a/Ln^ < 0) 
drive outward fluxes of main ions; inwardly peaked impurity profiles (a/Ln^ > 0) 
reinforce the inward fluxes of main ions (b). 



